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A UNIFORM AREA LAW FOR THE ENTANGLEMENT OF 
EIGENSTATES IN THE DISORDERED XY CHAIN 
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Abstract. We consider the isotropic or anisotropic XY spin chain in the presence of a 
transversal random magnetic field, with parameters given by random variables. It is shown 
that eigenfunction correlator localization of the corresponding effective one-particle Hamil¬ 
tonian implies a uniform area law bound in expectation for the bipartite entanglement 
entropy of all eigenstates of the XY chain, i.e. a form of many-body localization at all en¬ 
ergies. Here entanglement with respect to arbitrary connected subchains of the chain can 
be considered. Applications where the required eigenfunction correlator bounds are known 
include the isotropic XY chain in random field as well as the anisotropic chain in strong or 
strongly disordered random field. 
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1. Introduction 

1.1. Motivation: Towards Many-Body Localization. We are interested here in a bet¬ 
ter understanding of mathematical characterizations of many-body localization (MBL) in 
interacting quantum systems. This phenomenon has recently received strong attention in 
theoretical physics and quantum information theory, see e.g. [4, 5, 6, 10, 11, 17, 19, 22, 33, 
34, 36, 37, 38, 44, 45] and references therein. It is generally described as the absence of ther- 
malization or self-equilibration in a quantum many-body system, often due to the presence 
of disorder. For a detailed discussion of the current understanding of thermalization and 
many-body localization in the physics literature see also the recent survey [18], including 
the extensive list of references provided there. It is made clear in these works that there 
is no complete consensus yet on what physically constitutes thermalization and MBL, thus 
leaving a multitude of questions for further investigation and clarification. 

We will focus here on some of the criteria which by now are well accepted to be necessary 
characteristics of MBL and on studying these criteria for relatively simple models where they 
can be proven rigorously. Generally, and quite roughly, these criteria say that an interacting 
system, in suitable regimes such as weak interaction or large disorder, has properties similar 
to those of a non-interacting many-body system. The eigenstates of the latter are product 
states and thus have vanishing entanglement and spatial correlations. Also, its dynamics 
are trivial with no information propagating between particles. For an interacting system 
in the MBL phase one thus expects rapid decay of correlations and small entanglement 
of eigenstates, as well as absence of information transport (i.e. no or slow propagation of 
particle group waves). Note that this is a true many-body concept and to be distinguished 
from classical Anderson localization for a single particle, say, in a disordered environment, 

G. S. was supported in part by NSF grant DMS-1069320. 

1 



2 


H. ABDUL-RAHMAN AND G. STOLZ 


where ‘localization’ refers to the single particle configuration space. Significant differences 
between the many-body localized phase and single particle localization have been pointed 
out, e.g., in [4, 36]. 

Let us stress here that the term MBL should generally be reserved for properties which 
hold uniformly in the number of particles in the system (e.g. in the sense that relevant 
constants are bounded in the particle number). In this sense, many-body localization is to 
be distinguished from few-body localization, such as the known rigorous results for the N- 
particle Anderson model [13, 14, 2, 27], which do not yet allow uniform control in the number 
of electrons. In particular, results expected by physicists for the many-body Anderson model, 
such as MBL at low electron density or weak interaction strength, e.g. [5, 19], can not yet 
be shown rigorously. 

Rigorous mathematical results on localization properties of disordered many-body systems 
are so far essentially restricted to two models: Disordered harmonic oscillator systems [31, 32] 
and the XY spin chain in random held [28, 20, 41, 35]. These models are equivalent to free 
Boson systems and free Fermion systems, respectively, and thus can be studied in terms of an 
effective one-particle Hamiltonian. As a consequence, it is possible to deduce results on MBL 
from known localization properties of one-particle Hamiltonians such as the Anderson model. 
Of course, a long term goal must be to develop methods which allow to go beyond these 
simple models. In particular, an important challenge is to develop mathematical methods 
to study the disordered XXZ or Heisenberg chain, which reduce to the physically more 
interesting situation of interacting Fermion systems, see e.g. [4, 6, 33, 34, 37] for numerical 
results. Some progress towards understanding MBL for larger classes of spin systems in the 
presence or strong disorder has been made in [23], which, however, still needs to make an 
unproven assumption to limit the amount of level attraction in the system. Mathematically, 
describing the phenomenon of many-body localization and fully proving it for important 
classes of physical examples is a wide open held (even wider than for physics). 

Our more modest goal is here is to contribute to a more complete understanding of MBL 
properties of the two models indicated above, harmonic oscillators and the XY chain. 

For oscillator systems the reduction to an effective one-particle Hamiltonian is rather 
straightforward. In particular, the reduction works in any dimension and does not affect 
locality properties of the system. This has allowed to verify a quite complete list of MBL 
properties for disordered oscillator systems, including a zero-velocity Lieb-Robinson bound 
on information transport as well as exponential decay of correlations and area-law-type 
entanglement bounds for ground and thermal states [31, 32] (in each case requiring disorder 
averaging). 

The reduction of the XY model via the Jordan-Wigner transform to a free Fermion system 
is only possible for a one-dimensional chain of spins and, as an additional difficulty, introduces 
non-locality. Therefore rigorously known MBL properties for the XY model are more limited. 
A first contribution was made in [28] which established exponential decay of certain ground 
state correlations for the isotropic XY chain in random exterior held, using localization 
properties of the Green function of the effective Hamiltonian, in this case the Anderson 
model, proven via multiscale analysis. More general results have recently been obtained 
by Sims and Warzel [41], Using localization of eigenfunction correlators of the underlying 
one-particle Hamiltonian, these authors fold exponential decay of stationary as well as time- 
dependent correlations for larger classes of states, including general eigenstates as well as 
thermal states, for systems such as the XY chain which can be mapped to free Fermions. 
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Localization of eigenfunction correlators is also the key property of the effective Hamil¬ 
tonian used in [20] to show a zero-velocity Lieb-Robinson bound for the XY chain in random 
held, again after disorder averaging. An argument in [20], valid for a very general class of 
quantum spin systems shows that this implies exponential decay of ground state correlations 
up to a logarithmic correction in the size of the ground state gap (the more specific arguments 
used in [41] show that no such correction is required for the random XY chain). This in 
turn, by another general observation, leads to an area law for the ground state entanglement 
[7, 8]. Note, however, that [7, 8] consider deterministic systems and that some further anal¬ 
ysis would be required to show that the relation “exponential decay of correlations implies 
area law” carries over to the disorder averaged quantities. 

Here we prove an area law for the disorder averaged entanglement of all eigenstates of the 
XY chain in random field. In fact, we state this as a conditional result: If the corresponding 
effective one-particle Hamiltonian has localized eigenfunction correlators, then the rnany- 
body eigenstates satisfy an area law, uniform in energy (i.e. in the labeling of eigenstates), 
see Section 1.2 for exact statements. Applications of this result, see Section 1.3, include the 
isotropic XY chain but also some anisotropic cases. 

To our knowledge this is the first result on many-body localization for eigenstates of a 
disordered spin chain which goes beyond the ground state, the latter physically corresponding 
to the case of zero temperature. In fact, the uniform area law for all eigenstates found here 
corresponds to the infinite temperature limit where all eigenstates are equally probable. 
That the validity of an area law for all (or at least a large fraction of all) eigenstates is a 
necessary characteristic of the MBL phase was recently stressed by Bauer and Nayak in [6]. 
These authors provide a formal definition of the MBL phase, reflecting the idea of Fock space 
localization proposed in [5], and argue that an area law for the entanglement of eigenstates 
should generally follow as a consequence. It is natural to expect that localization at infinite 
temperature implies localization at finite temperature. However, our methods (in particular 
the use of the von Neumann entropy of reduced pure states to measure entanglement) do 
not directly extend to give useful entanglement bounds for mixed Gibbs states. 

Our result is closely related to recent work by Pastur and Slavin [35] who establish an 
area law for the ground state of disordered free Fermion systems. In fact, by mapping the 
XY chain to a free Fermion system via the Jordan-Wigner transform, we use the argument 
from [35] as one of the key ingredients to our proof. We extend this argument to general 
eigenstates and, to include the anisotropic case, also consider more general Fermion systems 
which do not conserve the particle number, leading to random block operators as effective 
Hamiltonians. It may also be worth noting that we consider entanglement with respect to 
arbitrary subchains Ai of the chain A = {1, ...,n} in Theorem 1.1 below (and not just 
left or right ends of A). This requires to work with a suitably chosen local version of the 
Jordan-Wigner transform, introduced in Section 3.3 below. 

Acknowledgement: It is our great pleasure to acknowledge many helpful discussions with 
Bruno Nachtergaele and Robert Sims on the contents of this work. G. S. would also like to 
thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, for support and 
hospitality during the program Periodic and Ergodic Spectral Problems in Spring 2015 where 
part of the work on this paper was done. 
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1.2. The Main Result. We consider the anisotropic XY spin chain in transversal magnetic 
field, given by the self-adjoint Hamiltonian 

n— 1 n 

(!) H ■= - Yl + ^) a f a f+i + ( X - <7j+i] - u i a f 

j =i j =i 

in "H = (^) ;gA C 2 , where A = {1,2,... ,n} for an arbitrary positive integer n. By af, crj 
and cr z we denote the standard Pauli matrices acting on the j-th component of the tensor 
product. 

The parameters [ij, 7 j and Uj describe the interaction strength, anisotropy and field 
strength, respectively, and we will think of them as the first n components of sequences 
of real-valued random variables indexed by j 6 N. While stronger assumptions will be 
needed in applications, for our main result we only will assume their uniform boundedness, 
i.e. that there exists C < 00 such that 


(2) sup \ + | 7 j| + \iSj\) < C almost surely. 

jm 

We will also assume that the disorder guarantees that 

(3) H almost surely has simple spectrum, 


in the sense that all its eigenvalues are non-degenerate. In the applications discussed in 
Section 1.3 below this will follow from an argument provided in Appendix A. 

As known since the work of Lieb, Schultz and Mattis [30], the XY chain Hamiltonian H 
can be mapped to a free Fermion system on H by the Jordan-Wigner transform. As a result, 
the explicit diagonalization of H reduces to the diagonalization of an effective one-particle 
Hamiltonian which can be written in form of the symmetric 2 x 2-block Jacobi matrix 


( 4 ) 


M := 


/ -u x a z fiiSi'yi) > 

-V2<J Z 

i^n—l^i, An—1) 

\ ^ n _i5 , (7 n _i) t -v n cr z J 


where cr z = 


1 0 

0 -1 


is the third Pauli matrix and 


( 5 ) 



For some more details on the reduction of H to M see Section 2 below. 

Our main result says that a suitable form of localization of the effective one-particle 
Hamiltonian M implies a uniform area law for the bipartite entanglement entropy of all 
eigenstates of the n-spin Hamiltonian H. To introduce the latter, consider an arbitrary 
subinterval Ai = {r,... ,r + £ — 1} of A and the bipartite decomposition "H = "Hi < 8 > ; H 2 , 
where 


Hi = 0 C 2 and n 2 = (g) C 2 . 

jeAi jgA\Ai 


(6) 
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For a pure state p in 7~L the entanglement entropy £(p) with respect to the subsystem Ai is 
given by the von Neumann entropy of the reduced state p\ = Tr Wa p, 


(7) 

£(p) = S(p i) := — Tr g(pi), 

where 


(8) 

, , f x logx, 0 < x < 1, 

9W = (o, x = 0. 


We choose log to denote the natural logarithm, as opposed to the dyadic logarithm used in 
this context in the information theory literature. The distinction is irrelevant for our work, 
as we will not keep track of universal constants. 

With E(-) denoting the disorder average, we will prove 

Theorem 1.1. Assume (2) and (3), as well as the existence of (3 > 2 and C < oo such that 

(9) E ll9(MWI ) 5 i + \T= W 

for all n G N and 1 < j, k < n. 

Then there exists C < oo such that 

(10) E ^sup£(p</,)^ < C 

for all n, r and l with l<r<r + £— 1 < n. In (10) the supremum is taken over all 
normalized eigenfunctions if of H and pp = |0)(0| is the orthogonal rank-one projection 
onto O0. 

In (9), g : M —> C may be any function such that |g(x)| < 1 uniformly, and g(M) is defined 
via the functional calculus for symmetric matrices. Using the block matrix notation from 
(4) above, the matrix elements g(M)jk are understood as 2 x 2-matrices. While any norm 
|| • || on the 2 x 2-matrices could be used in (9), we will work with the Euclidean matrix norm 
for the sake of definiteness. 

Bounds of the form (9) are generally referred to as localization of eigenfunction correla¬ 
tors. Their use in single particle localization theory originates from [1], see also [42] for an 
introductory survey. 

In particular, considering the functions gt(x) = e~ ltx , t 6 R, the bound (9) includes 
dynamical localization, uniform in time, for solutions of the Schrodinger equation if'(t) = 
—iMif(t). We note that requiring only power law decay for some (3 > 2 in (9) is a relatively 
weak form of eigencorrelator localization. In applications where (9) can be veriEed, one 
usually has much faster decay in \j — k\, see Section 1.3 below. 

The “area law” (10), giving an upper bound for the entanglement proportional to the 
surface area of the subsystem Ai (as opposed to the elementary bound in terms of its volume), 
is uniform not only in the size of the system A and subsystem Ai, but also applies uniformly 
to all eigenstates of H. It thus yields a form of many-body localization for H at arbitrary 
energy. 

Note that by a general result due to Hastings [21], an area law holds for the ground state of 
one-dimensional spin systems with a uniform (in volume) ground state energy gap. Not only 
is our work not limited to the ground state, but we also do not need any explicit assumption 
of the size of spectral gaps. Instead, we exploit single particle localization of the effective 
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Hamiltonian M as the mechanism to prove an area law. Of course, bounds on level statistics 
can be considered as being an implicit part of the underlying single particle theory. A price 
we have to pay for not having deterministic gaps between energy levels is the need to include 
disorder averages in (9) and (10). 


1.3. Applications. Here we discuss some concrete examples, i.e. assumptions on the ran¬ 
dom variables Pj , 7 j and Uj which guarantee (2), (3) and, most importantly, one-particle 
localization of M in the form (9). 

We will focus on the case of random magnetic field, where the corresponding one-particle 
localization properties are best understood, i.e. we choose pj — p £ M \ { 0 } and — 7 G ffi. 
to be constant, while the Uj are i.i.d. random variables with compactly supported distribution 
p. Thus ( 2 ) holds. Assuming, moreover, that p is absolutely continuous dp[y) = h(z/)dz/, it 
follows from Proposition A.l in Appendix A that (3) is satisfied as well. 

The localization property (9) is known for the following cases. 

(i) Let 7 = 0, so that 


( 11 ) 


H = 


-p 


B 


XX , Y Y \ 
CTj CF j+ 1 + (Tj a j+1 ) 


E 


•', n i 


is the isotropic XY chain in random field. Then all 2 x 2-matrix-entries of M are diagonal 
and thus M decomposes into A © (—A), where 


( 12 ) 


A 


[-"id \ 

p, 

\ H -u n J 


is the Anderson model on the finite interval A. If the density h is bounded and compactly 
supported, then it is known that, for some C < 00 and rj > 0, 


(13) E sup \g(A) jk \ < Ce-^~ k I, 

J 

see e.g. [42], This readily implies (9) with exponential decay in \j — k\. 

(ii) For the anisotropic case 7 ^ 0 , which does not reduce to the Anderson model, local¬ 
ization properties of one-particle Hamiltonians given by block Jacobi matrices of the form 
(4) have been studied more recently. 

A result by Elgart, Shamis and Sodin [16] covers the large disorder case: If the i.i.d. 
random variables j 6 N, have bounded compactly supported density, and if 1 7 = 
for A > 0 sufficiently large, then (9) holds with exponential decay in |j — k |. This is done 
in [16] through an adaptation of the fractional moments method to a class of random block 
operators which includes our model (4). 

If the magnetic field is strong enough to create a spectral gap for M around energy E = 0 
(uniform in the volume n and the disorder), then we can apply a result from [12]: Suppose 
that for some C > \p\ it holds that either > C for all j G N or Vj < —C for all j G N, 
resulting in the spectral gap (—[C — \p\),C — |/i|) for M. Then for every £ < 1 there is 
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C = C(£) and p = r/(£) > 0 such that 

(14) E [ sup ||g(M)j fc || J KCe-rt-^. 

\M<i / 

This is essentially what is shown in the proof of Theorem 7.2 of [12], while the results there 
are only stated in terms of dynamical localization for M, i.e. for the functions gt{x) = e~ ltx 
(but the argument covers general \g\ < 1). 

As discussed in [12], it remains an open problem if a bound such as (14) (or at least 
(9)) can be proven for M without the assumption of a zero-energy gap. A problem arises 
from the fact that the transfer matrix group (Furstenberg group) associated with M at 
E = 0 is not irreducible, so that many of the available dynamical systems tools for proving 
one-dimensional localization don’t apply. 

1.4. Notes on Content and Related Work. The rest of this paper is devoted to the proof 
of Theorem 1.1. As in most works on the XY chain, the starting point is the representation 
of the XY Hamiltonian 77, via the Jordan-Wigner transform, as a quadratic form in a set of 
Fermionic operators, see Section 2. A key tool is that eigenstates of the latter satisfy Wick’s 
rule (sometimes also referred to as being quasi free or gaussian ) and thus are fully determined 
by their correlation matrices with respect to the Jordan-Wigner Fermionic operators. 

More specifically, the crucial formula for determining the bipartite entanglement entropy 
of eigenstates is the trace identity (38) in Lemma 3.1 below. We will use it for the case 
where the state p in (38) is the reduction of an eigenstate of 77 to the subsystem 77 1 (which 
will ultimately be justified in Lemma 3.7). Thus the left hand side of (38) becomes the 
entanglement entropy of an eigenstate. We thus need to calculate the correlation matrix 
on the right hand side of (38), which will be accomplished in Lemma 3.2 and Theorem 3.4: 
Lemma 3.2 provides the central connection between the many-body Hamiltonian 77 and 
the effective one-particle Hamiltonian M, saying that the correlation matrices of eigenstates 
of 77 are given by spectral projections of M. Moreover, the proof of Theorem 3.4 shows 
that correlation matrices of reduced states are given by 27 x 27-restrictions of these spectral 
projections. 

Equipped with these tools we can then complete the proof of the area law (10) in Section 4. 
This is accomplished by using an extension of an argument from [35] to reduce the claim to 
the eigencorrelator localization bound (9). 

The trace identity (38) for Fermion systems seems to originate from [43], see also [29] for 
a more detailed presentation and [15] for a survey of subsequent results. In these works 
(38) was used as a tool in the study of the ground state entanglement for the XY chain in 
constant magnetic field, 

n— 1 n 

(15) - X)[(1 + 7)44k + (! - 7)44+1] - u Y. a r 

3 =1 3 =1 

In particular, the dimension reduction from 2 n to 2 n accomplished by (38) (or 2 £ to 27 for 
the subsystem) allowed numerical predictions for the limit of large subsystems Ai (and in 
the thermodynamic limit of an infinite chain A). These were later proven rigorously, see 
[24] for most complete results as well as references therein. Generally, it was found that the 
ground state entanglement for the anisotropic chain remains bounded in the size 7 of the 
subsystem, i.e. satisfies an area law. Phase transitions with divergent entanglement appear 


H. ABDUL-RAHMAN AND G. STOLZ 


for v = 2 and for the isotropic chain 7 = 0 . In particular, [25] shows that in the latter case 
the ground state entanglement grows as log l. 

Our work here differs in several significant respects from all previous works on entangle¬ 
ment in the XY chain and the related Fermion systems: We consider general eigenstates of 
the XY chain and not only ground states here, we work in finite volume (proving bounds 
which hold uniformly in the volume) rather than in the thermodynamic limit of an infinite 
chain, and, finally, use a modification of the Jordan-Wigner transform to be able to consider 
entanglement with respect to more general subsystems. Partly for these reasons, but also 
to make a number of tools from the physics literature accessible to a broader audience, we 
include a self-contained presentation of relevant parts of the theory of finite Fermion systems 
here. This includes thorough discussions of correlation matrices, Bogoliubov transforms and 
quasi free states. Some of this is contained in Section 3, with more background collected in 
Appendix B. 


2. Reduction of H to M via Jordan-Wigner 

The importance of the XY chain as a model in the theory of quantum spin systems goes 
back to the work of Lieb, Schultz and Mattis [30], where it was shown that the XY chain 
with constant coefficients (and initially without magnetic field) is an exactly solvable model. 
Their argument proceeds by using the Jordan-Wigner transform to reduce the XY chain to a 
free Fermion system. For the last half century this has turned the XY chain into a canonical 
toy model for quantum spin systems, frequently used as a first example to illustrate new 
concepts. 

It was understood that the methods of Lieb, Schultz and Mattis can be extended to include 
magnetic fields and to allow for variable coefficients /ij, 7 j and r]j. In the latter case the model 
is not exactly solvable, rather the diagonalization of H can be reduced to the diagonalization 
of the effective Hamiltonian M. 

Below we briefly summarize this reduction. A similar account with somewhat more detail 
(but different sign conventions) is provided in Section 3.1 of [20]. 

The local lowering and raising operators are 

(16) a := ^(u- Y - ia Y ) = ^ J q ^ , and a* = + ia Y ) = ^ , 

whose actions on the j'-th spin in Ft will be denoted by dj and a*. 

The Jordan- Wigner transform refers to the operators 

(17) ci := ai, cj := erf ... j = 2,..., n, 

which satisfy the canonical anti-commutation relations (CAR) 

(18) {cj,c* k } = 5 jk 1 , {cj,c k } = {c*,c* k } = 0 for all j,k = 1 , ...,n. 

Following the more general convention used in Appendix B below, we will refer to 

(19) C=(c 1 ,cJ,...,...,c n ,O t 

as the Jordan-Wigner Fermionic system. 

The XY chain H can be expressed in terms of the Jordan-Wigner operators as 

(20) H = CMC. 

Here the right hand side is interpreted in the sense of matrix multiplication of the col¬ 
umn C, the row C* = (c*,ci... , c*,c n ), and the scalar 2 n x 2n-matrix M given in (4). If 
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P is a permutation matrix which maps the canonical basis vectors ei,...,e 2 n of C 2n to 

('ll Cn+l; ^2; &n+2i ■ ■ ■ ■ &ni then 


( 21 ) 

where 


PMP t =: M = 


A B 
-B -A 



f-v i Hi 

\ 


/ 0 hHi 

\ 

(22) A = 

Hi 

Hn—l 

, B = 

~liHi '' • 

* ^/n— l/^n— 1 


V 


V 


^n—l^n —1 


We see A* = A t = A and B* = B l = —B, and thus M* = M l = M. 

'0 l' 


Transforming with the unitary 


1 0 


, we see that M is unitarily equivalent to —M, and 


thus, in particular, has spectrum symmetric to zero. It can be diagonalized by an orthogonal 
matrix W Let 


(23) 


W -=2 


v + u v-u 
v-u v + u 


where U and V are real orthogonal matrices associated with the singular value decomposition 
of S := A + B via 


(24) USV* = A := diag{Ai, A 2 ,..., A n }. 

Here 0 < Ai < A 2 < ... < A n are the singular values of S, i.e. the eigenvalues of (S*S)%, 
counted with multiplicity. Calculations show that 

(25) M = W‘(^ _° A ) W, WMW* = 0 (T 

where W = P 4 bFP is Bogoliubov, i.e. orthogonal and satisfies (92). Let 

(26) B := WC. 

By Lemma B.3 this is a Fermionic system and 

(27) H = CMC = B*WMW t B = B* 0 ( _° A . ) # 

3 =i ^ ' 

n n 

= E A AA - cc*) = 2 E - E ° 1 , 

3 =1 3 = 1 

where E 0 = y )™ =1 A j. Thus H has been written in the form of a free Fermion system. Let hi 
be the vacuum vector of the bj and 

(28) ip a = (bl) ai ... ( 6 *)" n H, a = (ati,ct 2 , ■•■,««) e {0, l} n 

the orthonormal basis of B associated with B , see Appendix B. The ip a form a complete set 
of eigenvectors for H with corresponding eigenvalues 2 Ylj-a=i Ai — -Ebj so that the spectrum 
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of H is 
(29) 


a(H) 


^2 Xj - ^2 Xj : a E {0,1} 


3. Entanglement Entropy of Eigenstates 

3.1. Correlation matrices and quasi free states. The goal of this section is to show 
that the entanglement of the eigenstates p a = |V’aXVul °f H with respect to the bipartite 
decomposition given by ( 6 ) can be expressed in terms of correlation matrices of size 

21 x 2£ and that these correlation matrices are restrictions of suitable spectral projections for 
the one-particle Hamiltonian M. This provides the crucial connection between properties of 
M and properties of H which will be exploited in Section 4 to prove Theorem 1.1. 

The expected value of an observable A E B{FL) in the mixed state p E B(FL) (i.e. p > 0 
with Tr p = 1) is given by 

(30) (A) p :=TrAp. 

Let D = (di, d*, d 2 , d * 2 ,..., d n , d * n )* be a general Fermonic system of Ft in the sense of 
Appendix B. The correlation matrix of the state p with respect to T> is defined to be the 
2 n x 2 n matrix 

(31) T v p := (VV*) p , 

with the row V* = (d\,d\ ... ,d*,d n ) and (VV*) p to be understood in the sense of taking 
expectations of each of the operator-valued entries of the 2 n x 2 n-matrix W*. 

Using the CAR, cyclicity of the trace, and trA* = trA, one checks that correlation 
matrices are generally of the form T = (r jk)\<j k <n with 2 x 2 -matrix-valued matrix elements 

(32) 

Here X and Y are n x n-matrices such that X* = X and Y* = —Y. Note that 

(33) Prp‘ = i( l2 „ + (; X F F)). 

In particular, correlation matrices, as defined here, are self-adjoint. We mention that in the 
physics literature frequently the Majorana Fermions dj + d*, —i(dj — d*), j = 1 ,,n, are 
used to define correlation matrices, which in this case become skew-adjoint. 

For / : A —> C define 

n 

(34) d(f) = Yji d i 

1=1 

and denote its adjoint by d*(f) = J2j fjdj- More generally, let 

(35) D(f,g) :=d(/)+cf(< 7 ). 
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We say that the state p is quasi free with respect to V if 


(36) 



r o, 

m 

< £(—1 ?{D t D,) p 

k=2 



if m is odd; 
if m is even. 


Here Dj is short for D(fj, gf) and a pair fj, gj : A —>■ C. Also note that here and throughout 
this paper operator products as on the left hand side of (36) are to be read from the right 
to left, i.e. n ™=i Dj = D m ... D 1 . 

For the case of even m, where (36) is a form of Wick’s rule, an iterative application shows 
that expectations D^ can be written as a sum of products of terms of the form 

(. D r D s ) p . We note that the resulting expression is known as the Pfaffian (denoted by pf) of 
the skew-adjoint m x m-matrix with entries Difr = ( D r D s ) p for 1 < s < r < n and 

appropriately extended by antisymmetry. The Pfaffian of an odd-dimensional skew-adjoint 
matrix is generally set to zero, so that (36) can be restated as 


(37) = pf W'"") ' 

As the expectations (. D r D s ) p are linear combinations of elements of the correlation matrix 
T®, we conclude that (^n”=i D^ is uniquely determined by the correlation matrix. By 

Lemma B.2 this determines (A) p for all A E £?("H), which in turn characterizes p. This 
yields the first claim of the following Lemma: 


Lemma 3.1. If a state p is quasi free with respect to V, then p is uniquely determined by 
the correlation matrix T®. Moreover, the von Neumann entropy of p is given by 

(38) S(p) = -Trplogp = -trrJ’logT^. 

Note here that Tr and tr denote the traces in B{fH) and C 2nx2n , respectively. The second 
claim is restated as Proposition B .8 in Appendix B, where a proof is provided. 


3.2. Correlation matrices of eigenstates. Next we calculate the correlation matrix of the 
eigenstates p a = \'ip a )('ip a \ °f H given by (28) with respect to the Jordan-Wigner Fermionic 
system C. 

Lemma 3.2. Assume that M has simple spectrum. Then, for each a G {0, l} n , the corre¬ 
lation matrix of p a with respect to the Jordan-Wigner Fermionic system C in (19) is given 
by 

(39) IL=Xa„(M), 

i.e. the spectral projection for M onto the set 

(40) A q := {A j : otj = 0} U {-A j : oq = 1}. 

Here { X 3 , —A j : j — 1, 2,..., n} are the eigenvalues of M . 
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Note that, in particular, p 0 = |f2)(f2| is the ground state projection for H and T c po = 
X{\i,...,x n }(M) = X(p,oo)(M). The difference between T^ q and T c pa is that A j is replaced by 
—A j for each site j in which a “particle is created” by b* in (28). 


Proof. The first step in the proof is to show that 


(41) 



0 


5 


where B is the Fermionic system (26). First note that, using the definition (28) and the 
anti-commutation properties of the bj, one has 


(42) 

b*i>a = j 

r o 

if a.j 

if aj 

as well as 




(43) 

bj'lpa = j 

r o 

( FpQ—ej 

if aj 
if dj 


Then statement (41) follows from the following calculations 


(44) Tr bjb* k p a = fy&fcV’a) = (frjV’o &/$>«) 

JO i f a 3 — 1 or a k = 1, 

\ (T"0a+ej j T'0 o ,_|_e fc ) 8j k if Qj OL k 0, 

And, 

(45) Tr bjb k p a = (^ a ,b j b k 'if a ) = {b*if a ,b k 'if a ) 

JO if o.j — 1 or d k = 0, 

\ ±{l/j a +e j , i>a-e k ) =0 if CKj = 0 and d k = 1, 


The general structure (32) of correlation matrices determines the remaining entries of T® 
and completes the proof of (41). 

We now prove (39). As C = W t B 1 by Lemma B.3, 


(46) 



= w t r B p w = w t Q) 

k=1 



0 

ba k , 1 


w. 


Now, since the eigenvalues of M are simple, i.e. — \ n < ... < —Ai < 0 < Ai < ... < A n , 
0£ =1 ^ a g’° j ^ ^ is the spectral projection for 0)) =1 ^ ^ onto A a . Therefore (39) 

follows from (25). □ 


Lemma 3.3. The eigenstates p a are quasi free with respect to the Jordan-Wigner Fermionic 
system C. 


Proof. Since B is a Bogoliubov transformation of C, by Lemma B.4(b) it suffices to prove 
that p a is Wick with respect to B. Let and = (cj) ai ... (d^) a "fl c be the vacuum and 
eigenstates associated with C. By Lemma B.4(a) we need to prove that pf' 1 := |) ffa' 1 | is 
quasi free with respect to C for all a. 
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This can be done by finding p'a explicitly. From the definition (17) of the c 3 we see that 


= 0 for all j. This gives that, np to a phase, D, c = 


and therefore 

(47) pf, C) = R>R| = 

Since for 1 < k < n we have 


is the vacuum for C, 




(48) 


n 


C j C j I C h 


VJ = 1 

it is easy to see that 
(49) 


3 = 1 

\ 


o o 

o l) _ ll c i c 3- 


n< 

3 = 1 


n 

3 = 1 

\ j 4 k / 


( 


c k c k c k c k 


C-jC 


3^j 

j = 1 

\ 3+ k ) 


C k C ki 


p? = n< 

3 = 1 




i!2c)Rin(%)“ j = n n 




3 = 1 


3 = 1 

Ct4 = 0 


j = 1 


i=i 


^a,.l 0 

0 <L, 


a ,,0 


Thus, as a product of diagonal states, R is quasi free with respect to C by Proposition B.6. 

□ 


By Lemma 3.1 we may conclude that Trp Q logp a = trT^logT^. This, however, is a 
trivial fact as p a is a pure state in 77 and Y c an orthogonal projection in C 2n and thus both 
traces vanish. The actual use of Lemmas 3.2 and 3.3 is that they provide a first step in 
establishing a similar trace identity for the reduced states (with non-vanishing entropy). 


3.3. Entanglement of reduced states. Restrictions of T^ q , which we have identified in 
Lemma 3.2 with spectral projections for M, can be used to express the entanglement of 
the states p a . Since p a is quasi free with respect to C we can apply the following general 
theorem for such states. Here, for any subinterval Ai = {r,..., r + £ — 1}, introduce local 
Jordan- Wigner operators {c^, j G on 771 as 

(50) R := a r , cf^ j = r + 1, ...,r + £ - 1. 

The operators 


(51) 


C i : = 



.,c 


(i) 

r+£— 1 ’ 


(c 


(1) 

r-\-l —1 


* 


t 


are a Fermionic system in 77i. 


Theorem 3.4. If p is quasi free with respect to C, then the entanglement entropy of p with 
respect to the bipartite decomposition 771 <g) U 2 is given by 

(52) £(P) = - tr T c p ] log Y c p ], 

where C\ is the local Jordan-Wigner Fermionic system on 77i. Moreover, 

(53) Y c p \ is the restriction ofY c p to span{e 2 j-i, e 2 j : j £ A^, 
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i.e. the 2£ x 2£-submatrix ofY c p consisting of the 2 x 2-matrix elements corresponding to 
The proof of this theorem fills the rest of this section. Let 

(54) Cj := tocf* ® j G Ai, 

be the extensions of the operators to TL. The c 3 are local observables for the subsystem 
Ai and differ from the Cj by products of ‘left-end’ Pauli matrices, 

(55) cj = of of ... of^Cj = Cj-o-fo-f ... of_,, j e Ai. 

To describe the difference between Cj and Cj further, for / : A x —* C we set 

(56) ~C r (f) = Y fjCj, C r (f) = Y Jj C V 

jeAi je Ai 

and denote their adjoints by c*(/) and c*(f). For pairs f,g : A x —>■ C we also write 

(57) C r (f,g) = c r (f) + c* r (g), C r (f,g) = Cr(f) + c*(g). 

From (55) and basic properties of the Pauli matrices, in particular (of... of _ x ) 2 = 1, we 
easily get 

Lemma 3.5. Let /, g, /i, f 2 , g 2 ■ A x ->■ C, then 

(a) C r {f,g) = of of .. .(Tr-iC r {f,g), 

(b) C r (f u gi )C r (f 2 , g 2 ) = C r (f u gi )C r (f 2 , g 2 ) 

With this we reach our first goal: 

Lemma 3.6. The correlation matrix Y c p \ of the reduced state p\ = Tr^ 2 p with respect to the 
local Jordan- Wigner Fermionic system C\ on TLi satisfies (53). 

Proof. That the upper right elements of the 2x2 matrix-elements in each correlation matrix 
in (53) coincide is seen as follows: For j, k G Ai, 

(58) (cfcfV = Tr (cf’ciV) = Tr (tfc' 1 * Tr Ms p) 

= Tr (Tr W2 c 3 c k p) = Tr ( CjC k p ) = ( CjC k ) p . 

Here we have used in the third step that 

(59) cf ] c^ Tt H2 p = Tr W2 (l ® cf c { p <g> 1) p = Tr Wz 5jC k p 

and in the fourth step that djC k = CjC k by Lemma 3.5(b). 

Identity of the other three elements of the 2x2 matrix-elements follows in the same way 
when replacing cp and/or cf by (c^f * and/or (cjf )*. □ 

In order to conclude that p\ can be determined from the correlation matrix Tf, we prove 
that is is quasi free with respect to C\. 

Lemma 3.7. If p is quasi free with respect to C, then the reduced state p\ = Tr^ 2 p is quasi 
free with respect to C\. 
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Proof. We need to prove (36) for pi-expectations of the operators 
(60) ) = c (1) (/) + (c (1) (p))*, c (1) (/) = jjcf 

je Ai 

on "Hi, where /, g : Ai —> C. The latter reduce to p-expectations of the operators C r (f,g) 
defined in (57): For any positive integer m and functions fj,gj : Ai —> C for 1 < j < m we 
have, 


(61) 


n cm (fps,) 

0=1 


■n n C (1) (fj, gi )Tr H2 p 


Pi 


0 = 1 
m 


mil 1.v= n«u) . 


0=1 


0=1 


which uses the general fact Tr^ 2 (i? ® I) A = B Tr^ 2 A for A e B(PL) and B e in the 

second step. 

Case 1: If m is even, then using Lemma 3.5 we get that 


(62) 


i =i i=i 


Note that the operators C r (fj,gj ) are of the form C(fj,gj ) where fj, gj are the extensions 
of fj, gj by zeros to A. Thus we can use that p is quasi free with respect to C, i.e. Lemma 
3.3: 


(63) 


n c(i) o.9j) 

0=i 


pi 


= pf (cy m) ), 

0=1 


where Cr p ' in ' > is the m x m anti-symmetric matrix with entries 
(64) [Cf* ra) ] jtk = (C r (f k ,g k )C r (f vgj )) p 

for 1 < j < k < n and properly extended by antisymmetry. Arguing as in (61) and (62), 
but now with only two factors, 


(65) 


(C r (fi,9i)Cr(fk,gt)) f = (C r (fj,g,)C r (h,3k)) P = 0C ,<1, (/i,®)C (1 >(A, 9t )>p 1 . 


Combined with (63) this shows that Wick’s rule is satished if m is even. 
Case 2: If m is odd, then by (61), (62) and Lemma 3.5 we get 


(66) 


n c<i) 0j.9j) 
0=1 


Pi 


m—1 

C r (fm,9m) JJ C r {fj,gj) 
3 = 1 
m—1 

= ( C r (fm, 9m) JJ C r (fj,gj) 

3 = 1 


m—1 


(uf a|... a^_ x ) C r (f m , g m ) JJ C r (fj, gf) 

3 =1 
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= (Vfa 2 z • • • af_i) PI Crifagj)^ ■ 

Now note that 

r— 1 

(67) a?<T% ... = H(c 3 + c*)( Cj - c*). 

3 = 1 

Thus erf erf ... erf_ 1 is an even product of operators C(fj,gj ) with suitable choices of /) and 
gj. Since p is quasi free with respect to C, we conclude that the right hand side of (66) is 
zero. □ 


Combining Lemma 3.7 and Lemma 3.1 gives 

(68) £(/>) = s(p 1 ) = -trr£;iogr6;, 

which completes the proof Theorem 3.4. 


4. An Area Law for the Eigenstates 


In this section we prove Theorem 1.1, our main result. 

Note first that, by the non-degeneracy assumption (3), it holds with probability one that 
all eigenprojectors pp of H are of the form p a = \i/} a )(i>a\- Therefore (10) is equivalent to 


(69) 


E 


ma x£(p a ) 

a 


< C < oo 


uniformly in n, r and £, meaning we can work with the Fermion basis given by (28). 

For fixed a, set Fj := T^^ the correlation matrix of (p a )i '■= Tr p a with respect to C\. 
r! is the restriction (53) of T := T c pa = XA a (M). Let the eigenvalues of Ti be 1 — 
j — 1 ,... ,£, with 0 < < 1/2. 

The following is a calculation essentially taken from [35], extended to the more general 
type of correlation matrices needed here to cover quadratic Fermion Hamiltonians of the 
form (20). 

By Theorem 3.4 we have 


(70) 


£(pa) = 


< 


- tr T 1 log Ti 

e 

~ 51 (& lo S& + ( x “ &) - £j)) 

3 =1 

e 

2 log2 

1=1 

i og 2tr(r 1 (i —r i)) 1/2 , 


where we have used the elementary inequality 

(71) — xlogx — (1 — x) log(l — x) < 21og2A/x(l — x), for 0 < x < 1. 
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The Peierls-Bogoliubov inequality, see Section 8.3 in [40], says that 

m 

(72) tif(A)>J2f(A sj ) 

3= 1 

for any convex function / and m x m hermitian matrix A. Using this with f{x) = —yfx as 
well as the elementary inequality yfx + yjy < \f2yjx + y, we may further bound (70) by 

(73) S{p a ) < V2\og2 Y, (tr(r!(l -T,))^ 2 , 

jeAi 

where matrix elements should be understood as 2 x 2-matrices. 

Now, since T is an orthogonal projection, we use T 2 = T with block matrix multiplication 
to get, 

(74) Tjj = Fjj + TjkTkj + r jkTkj- 

k e Ai fceA\Ai 

3 + k 

Then, for j £ A 1; 

(75) (Ti(l — Ti))^- = Tjj(l — Tjj) — TjkTkj = TjkTkj— Tjk(Tjk) , 

k e Ai fceA\Ai fceA\Ai 

fc + 3 


where symmetry of T was used (note that Ti is real-valued, see (46), where W is real). 
Inserting this into (73), using yjx + y < yfx + y/y as well as tr Tjk(TjkY < 2||r ?fc || 2 , we find 

(76) %)< 2 log2 ^ ]T |ILt||. 

jeAl keA\Ai 


Maximizing over a and averaging gives 



Now, by assumption (9), 


(78) 


E sup 


[XA a (Af)] 



< E 




< 


C 

1 + 1 j ~ k\P 


for some (3 > 2. We have 


( 79 ) 


1 -j- \j — k\P 

ieAi fceA\Ai 1 


' 1 + \j — k\P 

£ oo ^ 

2 Y1 i + (fc_ 

j =i k=e+1 y 

„ ^ i i 

£ Ukiu Vl + 2((’-j)'> W + 2(i - W 
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< 



VT+W 


< oo. 


This gives the uniform boundedness of (77) in n, r and £, and thus completes the proof of 
Theorem 1.1. 


Appendix A. Non-Degeneracy of Eigenvalues 

Here we will prove that the non-degeneracy assumption (3) holds under the condition that 
the random variables u y j £ N, are i.i.d. with absolutely continuous distribution and that 
Hj, t j, j £ N, are independent of the In particular, this covers the applications discussed 
in Section 1.3. In this case, almost sure non-degeneracy of the eigenvalues of H will follow 
from Proposition A.l below. 

Towards this, let A is a hermitian n x n-matrix, B an anti-hermitian n x n-matrix, A(i>) = 
A + diag(n) for v — (zq, ..., i/ n ) £ R n and 

< 8 °) -!(,))■ 

As M(v) and —M{v) are unitarily equivalent, M(v) has n pairs ±A j, j = 1, ...,n, of 
eigenvalues, where we may choose 0 < Ai < ... < A n . As discussed in Section 2, the 2 n 
eigenvalues of the corresponding many-body Hamiltonian H(v) = C*P t M{y)PC are 

(81) E a := Y. ^ - E A 1 

j:aj=l j-Oij =0 

for multi-indices a = (oq,..., a n ) £ {0, l} n . 

Proposition A.l. For Lebesgue-almost every v = (zq,..., u n ) £ M" the 2 n eigenvalues 
{E a : ol £ {0,1}"} of H{v) are pairwise distinct. 

Remark A.2. Pairwise distinctness of the many-body eigenvalues E a readily implies pairwise 
distinctness of the one-body eigenvalues, i.e. Proposition A.l yields that 0 < Ai < ... < X n 
and, in particular, invertibility of M(v) for Lebesgue-almost every v £ M n . 

The elementary proof of Proposition A.l proceeds in two steps: First we show the existence 
of one v' £ M n such that the 2 n eigenvalues of H(u') are pairwise distinct. Then we use 
analytic perturbation theory to show that this property extends to H{y) for Lebesgue-a.e. 
v £ M n . 


Step 1: Consider the functions fp(x) := YTj=\ ftj x j on P e {—1, l} 7 ^. There exists 
a nullset N C W 1 such that the numbers fp(x), /3 £ { —1, l} n , are pairwise distinct for 
all x £ M n \ N (first observe that for each pair /3 ^ (3 there is a nullset such that 
fp{x) 7^ fp{x) for all x £ M n \ , then take the union of these nullsets). 

Fix v" £ M n \ N with 0 < u'f < ... < u” and let 5 := min^ \fp(v") — fp{y")\. Let 


D : = 


(A 

B\ 


-a) 


and z/ := AnD+1 p ", S o that 


min | fp(v') - fp(i/) | = 4 nD + 1. 


(82) 
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Suitable choices of /3, /? show that \v' k \ > 2D and \u' k — v' e \ > 2D for all k i. Treating 
M{y') as a perturbation of diag(VJ,..., u' n , —is [,..., —u' n ), we see that each of the 2 n intervals 

(83) [uj — D, z/ + D], [—z/ — D, —Uj + D], j = l,..., n, 

contains exactly one eigenvalue of M{y'). Denote these eigenvalues by ±A'-, j = 1 ,,n. 
For the corresponding eigenvalues E'p := fp( A' 1; ..., X' n ), (3 e {1, — l} n , of if (z/) we get 

n 

Dft - W' 

3 = 1 
n 

_ &)W - v 'i) + 

3 = 1 

n 

> \UW)-W)\-^\K~^\ 

3 =1 

> |//3(^') — /^(^')l — 2rziD > 0 

if /3 /3. Thus if (z/) has non-degenerate spectrum. 

Step 2 : We now apply analytic perturbation theory iteratively to each of the parameters 
I'D • • •, v n in if( v). 

Fix z/ 2 , ..., z/, the last n — 1 components of z/ found above. Then 

(85) ... X) = M + i/i-Bi 

for selfadjoint Ai and f?i. Thus by Theorem 11.6.1 of [26] there are real analytic functions fk, 
k — 1,..., 2 n , such that {fk{y i) ■ k — 1,..., 2 n } are the eigenvalues of H(u i, v ' 2 ,..., i/) for 
each it e R. As H( z/) has non-degenerate spectrum, the numbers ) : A; = 1,..., 2 n j 

are pairwise distinct. Analyticity of the functions fk implies the existence of a nullset A^ 1 ) C 
R (in fact a countable set) such that the numbers {/&(it) : = l,...,2 n } are pairwise 

distinct for each ^ gK \ Ad 1 ) (each pair of functions fk and f k , k =fk, coincides at no more 
than countable many points). 

For any such zr, we can now use analyticity in v 2 (given that if (it, v 2 , *4, • • •, v' n ) = 
A 2 + z/ 2 f> 2 ) to get the existence of a nullset Ah 2 )(zr) in R such that if (zr, ^2, ■ ■ ■, z^) has 

pairwise distinct eigenvalues for all zr e R \ Ad 2 )(zr). ft follows that the eigenvalues of 
H (zr, zr,1/3 ,..., z/) are non-degenerate for all zr, zr with ^ G R \ A^ 1 ), f 2 e R \ N^ 2 fui), 
and thus for Lebesgue-a.e. (zr, v 2 ) € R 2 . 

Iteration of this argument leads to Proposition A.l. 

Appendix B. Finite Fermionic Systems and Bogoliubov Transformations 

In this appendix we collect some further background from the theory of finite Fermionic 
systems, which was used in Section 3 above. Most of this is well known in theoretical physics 
and can also be found in the mathematical physics literature, e.g. [9] or [3]. We include a 
self-contained presentation of this material here, partly for the convenience of the reader, 
but also to ensure that these tools are available in the generality required here. 

In a Hilbert space T~L of dimension dim T-L = 2 n , we call 

(86) V = (d u dl,d 2 ,d’ 2 ,.. ,,d n ,d’ n )‘ 


(84) 


iy> - 
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a Fermionic system if the operators dj G B(fH) and their adjoints satisfy the canonical 
commutation relations (CAR) 

(87) {dj,d* k } = 6 jk 1, {dj,d k } = {d*,d* k } = 0 for all j, k = 1,..., n. 

The intersection of the kernels of the dj is one-dimensional, i.e. they contain an essentially 
unique normalized vector 0, referred to as the vacuum vector, from which an orthonormal 
basis of Li is found as 


(88) <ft a =* K) Q1 ... «) an O, ae{ 0 ,ir, 

see for example [39]. 

It is easy to see that the operators {d\, d \,..., d n , d * n } in a Fermionic system are linearly 
independent (consider a general linear combination of these operators and calculate its anti¬ 
commutators with all dj and d*). Next we state two other basic properties of finite Fermionic 
systems. 


Lemma B.l. Let Li and Li be 2 n -dimensional Hilbert spaces and V = (d±,d\,... , d n , d*)* 
be a Fermionic system in Li. Then 

(89) V = (d 1 ,dl...,d n ,d* n ) t 

is a Fermionic system in Li if and only if there exists a unitary operator U : Li Li such 
that U*djU = dj for all 1 < j < n. The unitary operator U is characterized by U(j) a — cj) a 
for all a G {0, l} n , where {ft a } and {<ft a } are the ONBs associated with V and V through 
( 88 ). 


Proof. For U as above and any j we have 


(90) d*U(ft a 



0 , if atj 

(_l)Ei=i«fc0 a+ei = (-l)Efc = i^[/0 Q+e ^ if a . 


Thus d* = U*d*U and U*djU = dj for all 1 < j < n. 
The converse is straightforward. 


i; 

o. 


□ 


Lemma B.2. Let V be a Fermionic system in Li and A be the -k-algebra generated by the 
components ofV. Then A = B(Li). 

Proof. With <p a as in (88) it suffices to show that for every pair a, (3 G {0, l} n there exists 
an operator A a ^ G A, such that A a ^(f) a = (ftp and A a ^(fta = 0 for a ^ a. This operator is 
explicitly given by 

(91) := III (‘C-W') 111 "" ) . 

We omit the somewhat tedious calculations needed to verify this. □ 

A matrix W G C 2 " x2n is called a Bogoliubov matrix if W is unitary and 

(92) WJW t = J, where J = (a x )® n 

The reason for using this terminology is that, for the finite Fermionic systems considered 
here, Bogoliubov matrices implement Bogoliubov transformations: 
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Lemma B.3. Let V be a Fermionic system in FL and W € C 2nx2n . Then 

(93) V := WV 

is a Fermionic system in FL if and only if W is a Bogoliubov matrix. In this case, the 
correlation matrices, as defined in (31) above, are related by 

(94) rf = WT v p W* 
for all states p in FL. 


Proof. Note first that (94) follows from a simple linearity argument. Let W be a Bogoliubov 
matrix and T> a Fermionic system. The latter means 

(95) VV* + J (' DV*y J= l 2n . 

Note that, given unitarity, the condition (92) is equivalent to JW = WJ and W t J = JW*. 
Thus 


(96) 


VV* + J(VV*yj = W (VV* + J(VV*yj) W' = W 1 2 „ W* = 1 2 „, 


so T is a Fermionic system. 

By performing a simple change of basis, the converse can be restated as: If T> and T> are 
two Fermionic systems related by 


0 I 

i or 


( 97 ) [d u d 2 ,..., d n , (Jo*,..., {d n y J = w (di, d 2 ,..., d n , (di)*, ..., (d n yy 

then W is unitary and satishes 

(98) W.J X W l = ./,, where ./, : = 

To show this, write W in block form, 

(99) W = 

From the transformation (97) we get that for j = 1,2,..., n, 

n n 

(100) dj = y^(iFjfcdfc + Lj k d * k ), d* = ^( Mj k d k + Nj k d* k ). 


K L 
M N I ' 


k=1 


k=1 


By taking the adjoint of the left hand side of (100) and comparing it with the right hand 
side, using that the elements of a Fermionic system are linearly independent, we get 

(101) N jk = K jk , M jk = L jk , for all j, k = 1, 2,..., n, 
meaning that 

(102) M = L and N = K. 

Next, we will prove that W is a unitary matrix, that is the rows are an ONB of C 2n . Using 
(100) and (102), a calculation starting with {dj, d*} = 1 proves that the rows of W are unit 
vectors. Next, for j ^ k , two calculations starting from {dj,d* k } = 0 and {dj,d k } = 0 show 
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that any two rows of W are orthogonal. Thus W is a unitary matrix, which means that 


1 . In terms of the blocks of W this means that 


WW* 


KK* + LL* = 1, KL l + LK l = 0, 
LI<* + KL* = 0, 'LL 1 + KK* = 1. 


(103) 


One checks that this is equivalent to (98). 


□ 


In the situation above we say that T> is a Bogoliubov transformation of T>. One can easily 
see that being related by a Bogoliubov transformation is an equivalence relation between 
Fermionic systems. 

As in Section 3.1 we say that a state p on LL is quasi free with respect to a Fermionic 
system V if expectations of products of operators D(fj,gj ) satisfy (36). 

Lemma B.4. Let T> and T> be Fermionic systems on Ft and assume that p is quasi free with 
respect to V. 

(a) Let U be the unitary relating V and V as in Lemma B.l. Then UpU* is quasi free 
with respect to T>. 

(b) If V is a Bogoliubov transformation ofV, then p is quasi free with respect to V. 
Proof, (a) For any positive integer m, 



(104) 


From this it is straightforward to check that UpU* is quasi free with respect to T>. 


(b) Since the dj and d* are linear combinations of the dj and d*, one can see that, for 
given f,g : {1 ,... ,n} —► C, 


(105) 


D(f,g) = D{h,r) 


for some h, r : {1, 2,..., n} —> C. With this the claim follows easily. 


□ 


The following result states that, up to unitary equivalence, quasi free states are determined 
by their correlation matrices. 

Lemma B.5. Let p and p be states in 2 n -dimensional Hilbert spaces H and PL, respectively, 
which are quasi free with respect to the Fermionic systems T> and T>, respectively. If 


(106) 

then p and p are unitary equivalent. 



Proof. Using notations as in Section 3.1 (for both V and V), equality of the correlation 
matrices gives the equality of the correlations 


(107) 


(D,D k ) p = (D,D k )p for all 1 < j < k < n. 


Thus for any positive integer m, 
(108) 




As p and p are quasi free with respect to V and T>, this implies 
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u d > 

' UpU* 

where we used Lemma B.l which provides a unitary U such that dj = U*djU for all 1 <j< 
n. By Lemma B.2 this implies (A) p = (A)u p u* and thus 

(110) p = UpU* 

□ 

Our final goal in this Appendix will be to prove the trace identity (38). We will first 
prove this identity for diagonal product states with respect to the Jordan-Wigner Fermionic 
operators C given in (19). The general result (38) will then be reduced to this special case. 

Proposition B. 6. Let p( dia d) g #(0" C 2 ) be the product state given by 


( 111 ) 


P 


( diag) _ 


3 = 1 


Vi 0 

0 1 - TjjJ ’ 


where rjj E [0,1] for 1 < j < n. Then p( dm 9 ) j s quasi free with respect to the Jordan-Wigner 
Fermionic system C. Also, 


( 112 ) 

and 

(113) 


rc 

x p (diag) 


© 

3 = 1 


1 - Vj 0 

0 rjj 


Tr p ( dia 9) log = tr T c , diag) logP 


pyu,tu,y j o p(diag) 


The strategy of the following proof is essentially the ‘classical’ argument for the proof of 
Wick’s rule for thermal states in free Fermion systems, e.g. [9]. 

Proof. We first note that the formulas (112) and (113) follow from explicit calculations. For 
the rest of the proof we will drop the superscript (diag). For odd values of m in formula 
(36), the result follows from proving 


(114) 


n 

W-I 




= o, 


where the symbols ff rj stand for ff or nothing. We will assume that the product of c^’s is 
not zero, otherwise there is nothing to prove. Since 

(115) aa*a = a, a*aa* = a*, and a*a z = ±a # 

it is easy to see that 


(116) 


=±®a, 

3 = 1 j =1 


where Aj e { aj , a*, aja*, a*aj, a z , 1}, and since m is odd there exists j 0 E {1,2,..., rz.} such 
that Aj 0 E {aj 0 ,a* o }. Then 


(117) 


Tr iL ; * r v=±n ( tr ' 4 

3=1 3= 1 


rjj 0 

0 1 -rjj 
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which vanishes because 
(118) 


tr a 


# 


, ° i=o. 

0 1 - Ik, 


The proof for even m is more involved. First, assume that rjj ^ {0,1} for all 1 < j < n. 
Note that 


(119) 
because 

( 120 ) 

( 121 ) 

Also, 

( 122 ) 

where 

(123) 


Ckp 


Vk 


~ pC-k i 


a 


1 - Vk 

Vj 0 ^ = (Vj 0 \ = (Vj 0 

0 l~ Vj ) Vo -(1 -Vj)) V° 1 ~ r ) 

Vk 0 A = f 0 o\ _ Vk hk o 

o 1 - Vk) \Vk oj 1 - rj k V 0 1 - Vk 

n n 

c (f)p = EAc fc p = pX> i Vk Ck = P c ( D d)’ 

t — Vk 




k =1 


fe=l 


:= diag {£,■ : j = 1, 2 ,..., n} , ^ := 


Vj 


1-Vj 


c*{g)p = pc*{D^ g ) 


Similarly, 

(124) 

The CAR imply 

(125) (c(/), c(g)} = 0, {c(/), c*(g)} = (/, g)p. 

With Cj := C(fj , gj), j = 1,..., n, write 


(126) 


n a 

\j =i 


n g ) 

Kj=2 


P ' P 

For the first term we calculate, using (122) and cyclicity, 


n<Rc*( S i 

0=2 


( i27 > ((n c d c (/ 

\j=2 


1 ) J = 1 c 

P 


Wi)U c i 


m— 1 


= ( Ci 


(^/i)c m n c '^- 


1=2 


P 

m—1 


J=2 


P 

m—1 


ic(Dj ,),c m } n g > - (c m c(B e /i) n g 


1=2 


1=2 


P J ' P 

Im— 1 \ / m—1 \ 

= (SAhlf U c i - (c„c(c { /0 nCj) . 

' 3-2 / „ \ 3-2 / „ 

where also (125) was used. We proceed by applying this argument iteratively to the second 
term in (127), commuting c(D^f) with each of the Cj, and eventually conclude that (127) 






AREA LAW FOR THE DISORDERED XY CHAIN 


25 


coincides with 


(128) 


EMm/biftth n o) -(mo c(d £ /i 


k =2 


3 = 2 


0=2 


Defining fi := (1 +Dg)/i, the outcome of this calculation can be rewritten as 


(129) 


IO cCfO) =£(- 1 )‘< D «( 1+0 «)" I /i’9*W J] O 


0=1 


k=2 


3 = 2 
3 ¥= k 


Now 

(130) (D ? (l +-Df) _1 /i,0fc)<2 = (D v fi, gk), where D v := diag{r/j, j = 1,2,.. .,n} 
and, as we can explicitly check that (c*Cj) p = r)j, the latter coincides with 

(131) (cW/i)) = {c*(9k)c(fi)) p 


= (D r J\,g k ). 


3 = 1 


Thus equation (129) becomes 


(132) 


iq 4) =^(-!)‘(c t c(/,)>,( n O 


0=1 


k=2 


3 = 2 

J 7 k k 


By applying similar steps (which we omit) to the second term of (126), and introducing 
()\ = (1 +D^ 1 )gi in the process, we get 


(133) 


Ioh*®,)) =E(- 1 )‘<Dc-(ji))h n o 


0=1 


k=2 


3 = 2 
3 + k 


By substituting the results (132) and (133) into (126), we obtain Wick’s rule 

/ m \ m / m 

(134) (no) =E(- i )‘<c t c 1 ) P ( n o 


0 = 1 / p 3 = 2 


i ••• 2 

3 + k 


for the case of even m (after renaming fj and g 3 as fj and gf). 

Finally, in the general case where rjj G [0,1] for all j — 1, 2,..., n, there exists a sequence 


(135) 


Pn — 


(n) n 

m 0 


ML V 0 '-'Z 


(n) 


where ^ {0,1} —>• 77^ as n —>■ 00 and thus p n —> p. Now the fact that p n is quasi free 


with respect to C carries over to the limit. 


□ 
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We can now prove the following fundamental relation, which was stated earlier as the 
second part of Lemma 3.1. 


Lemma B.7. LetT> be a Fermionic system and p is self adjoint onB{FL), then the correlation 
matrix Y c p has a symmetric spectrum around | and it is diagonalizable by a Bogoliubov Matrix. 

Proof. It is enough to show that 
(136) 


T := 2r® - 1 


1 1 


-i i 


is diagonalizable by a Bogoliubov matrix. Let us define the unitary matrix Q := 

One can check easily that —iYlTYl* L is real anti-symmetric matrix, i.e. I = — r*. From 
the spectral theory of anti-symmetric matrices, there exists an orthogonal matrix O such 
that 


(137) 


r = o® 

3 = 1 


0 A, 
-A, 0 


0\ 


where A,- > 0 for j = 1,2 ,... ,n. By diagonalizing we get 


(138) 

Thus 

(139) 


QTYl* =iT = -O 
2 


r=W,| 0 -A, 


l —l 

1 1 


©n n 


© 

3 =1 


Aj 


0 


0 -Xj 


—i 1 
i 1 


Of 


W* where W = -=FL*0 

72 


—i 
1 


Note that W is unitary and one can easily check that (92) is satisfied. This proves that W 
is a Bogoliubov matrix. Then, (136) and (139) imply that 




!+Aj 

2 

0 


0 

1—Aw 


(140) 

j-i v u ^ 
which shows that the spectrum is symmetric around A. 




□ 


Proposition B.8. Let V be a Fermionic system in Ft and p a state in Ft which is quasi free 
with respect to T>. Then 


(141) 


Tr plog p — Ir Iy log 


Proof. Lemma B.7 implies that there exists a Bogoliubov matrix W and real numbers r/j 
such that 


(142) 




= : D 


V 


In fact we have q 3 G [0,1]: For the Fermionic system V = WT>, we get by (94) that = D v 
and thus rjj = Tr d*d 3 p G [0,1] because d*d 3 is an orthogonal projection. 

Consider the product state p( dm9 '> of Proposition B.6 with this choice of the numbers rp. 
An explicit calculation shows T c p(diag) = D v , which is therefore equal to Also, p < ' dm F 

is quasi free with respect to C by Proposition B.6 and p is quasi free with respect to V 
by assumption and Lemma B.4(b). Thus, using Lemma B.5, we get that p and are 
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unitarily equivalent. As we have also shown that T® is unitarily equivalent to T c p(diag) , (141) 
now follows from (113). □ 
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